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1.0 SUMMARY 


An integrate autopilot/autothrottle has been developed for implementation on the NASA TSRV B737 
airplane. The system was designed using a Total energy concept’. Certain maneuvers require a net 
energy change (e.g., an increase in velocity or height) which can be obtained using the throttle. In gusty 
conditions, the system may only require a redistribution of energy obtained from using the elevator. 

The system is intended to achieve: 

(1) Fuel efficiency by minimizing throttle activity. 

(2) Low development and implementation costs by designing the control modes around a fixed 
inner loop design. 

(3) Maximum safety by preventing stall and engine overboost. 

The control law was designed initially using linear analysis techniques and developed using a simplified 
nonlinear simulation. The system satisfied the design requirements for stability and cross coupling errors 
during maneuvers. The wind shear criterion ( < 5 kt vel error) was satisfied during approach while 
throttle activity was minimized in cruise. 

The control law was developed to include additional control modes (e.g., glide slope capture) and modelled 
on a nonlinear aircraft simulator (i.e., Harris simulator) to confirm the original design results and 
validate the system for additional flight conditions. 
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2.0 INTRODUCTION 


A new airplane control system must demonstrate more than just good dynamic response to be seriously 
considered as a replacement for existing systems. In the last decade the dramatic increase in the price of 
oil, cost of manpower, and amount of air traffic have meant additional constraints on developing a new 
system. The result of these constraints is that a new system must be designed for: 

(1) Fuel Efficiency 

(2) Minimum Development, Implementation and Certification Costs 

(3) Maximum Safety 

The integrated autopilot/autothrottle has been designed with these constraints in mind. To maximize 
fuel efficiency the integrated autopilot/autothrottle has been designed using the concept of a 'total energy 
control*. 

Consider an airplane flying at some datum condition with an energy level consisting of the sum of the 
kinetic energy (KE) and potential energy (PE), where the KE is associated with speed and the PE with 
altitude. In such a condition certain maneuvers may require a net energy change (e.g., an increase in 
velocity or height) which can be achieved by control of thrust using the throttle. In gusty conditions or 
during the execution of a double maneuver (i.e., decrease in height coupled with increase in velocity) the 
net energy change may be small, but the system will require a redistribution of energy which the elevator 
provides via redistributing energy by trading KE for PE or vice versa. 

This concept produces a 'coordinated* response (i.e., the throttle and elevator working in unison to the 
command inputs or disturbance), whereas in conventional autopilots and autothrottles (i.e., in which the 
elevator controls height and the throttle speed) the system can produce an 'uncoordinated* response 
requiring excessive throttle motion. 

This system has been designed with a constraint of minimum complexity to lower costs, and thereby 
reduce software development and verification compared with a conventional system. In achieving this 
aim, the system has been developed around a generalized inner loop fixed for all control modes, and an 
outer loop reconfigured for each required control mode. Furthermore, to maximize safety, the engine 
controller has been designed to include EPR limit protection, while the system incorporates angle of 
attack limit protection. 
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This report documents the development of the integrated autopilot/autothrottle and includes: 


( 1 ) Linear Design of the System 

(2) Development of Additional Features to Improve Performance in Nonlinear and Limiting 
Conditions 

(3) Design of the Engine Controller 

(4) Performance Evaluation of the System Using the Harris Nonlinear Flight Simulator. 

Detailed system diagrams are included herein to facilitate system software specifications and future 
flight test of the system. 
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3.0 NOMENCLATURE 
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Longitudinal acceleration (through center of gravity) 
Elevator command 
Elevator deflection 

Commanded Engine Pressure Ratio (EPR) 

Minimum EPR 

Total Energy Controller Error 

Total Energy Controller 

Height 

Height rate 

Height acceleration 

Barometric height 

Height command 

Inertial height rate 

Engine controller forward path gain 

Engine controller feedback gain 

Velocity error gain 

Height error gain 

Error gain 

Current pressure, pressure at sea level 
Longitudinal gust velocity 
Flight path acceleration 
Air referenced velocity rate 
Calibrated Air Speed (CAS) 

Equivalent Air Speed (EAS) 

Maximum EAS 

Velocity rate/acceleration error 
Thrust derived velocity rate 
Inertially referenced velocity rate 
True Airspeed (TAS) 

Throttle command 

Throttle level angle 

Weight, weight at sea level 

Vertical Gust Velocity 

Angle of attack 

Flight path angle 

Flight path angle error 

Elevator deflection 

Change in Thrust 

Airplane body angle, body rate 

Complementary filter (V I/AIR ) time constant 

Complementary filter bandwidth 

Natural frequency 

Damping ratio 

Thrust change 

Flight path angle change 
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4.0 DESIGN AND PERFORMANCE REQUIREMENTS 


The following design and performance requirements were imposed on the design to assure satisfactory 
performance of the final control law: 

(1) All vertical control modes should be integrated into a single control law concept. The 
system should be designed around a fixed inner loop configuration with outer loop control 
modes generating control signals that are compatible with the inner loop command 
inputs. 

(2) Cross coupling control errors should be minimized and the system should give a balanced 
performance over all flight conditions. The maximum velocity error should be less than 
1 kn for path control, while the maximum height error should be less than 20 ft for speed 
control. 

(3) There should be no overshoot of response for step command inputs. The damping ratio of 
dominant poles should be greater than 0.7. 

(4) Tight speed control is required in conditions of wind shear (max speed error < 5 kn). This 
criterion is considered vital at low speed and altitude conditions. 

(5) Throttle activity should be minimized particularly in cruise conditions. However, during 
approach when wind shear becomes a dominant consideration, the requirement to 
minimize throttle activity becomes a secondary consideration. 

(6) The design should incorporate safeguards to prevent stall and overboost of the engine. 
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5.0 LINEAR DESIGN 


5.1 INTRODUCTION 

The concept of 'total energy control* evolved during the evaluation of two preliminary control law 
schemes. The basic configurations and reasons for rejection of these configurations are presented here. 

5.1.1 Configuration 1 

Configuration 1 was concerned with integral control of speed through the throttle with proportional and 
rate control of speed plus integral control of altitude through the elevator. This system was found to 
require substantial modification to achieve system stability. In addition, performance characteristics 
were found not to be good. 

5.1.2 Configuration 2 

Configuration 2 was concerned with integral control of speed through the throttle and integral control of 
altitude through the elevator with proportional and rate control of speed and altitude through the throttle 
and elevator. After some modification from the original concept, this design yielded good performance in 
turbulence, despite that the cross coupling and windshear performance were poor. 

Nevertheless, Configurations 1 and 2 did not lend themselves to the single basic control law concept and 
required considerable software to provide all the additional control modes needed in a practical autopilot 
design. Therefore, a third control configuration was developed. 

5.1.3 Configuration 3 

Configuration 3 was concerned with integral control of altitude through the throttle and integral control 
of speed through the elevator with proportional and rate control of speed and altitude through the throttle 
and elevator. 

This third concept was similar to the technique used by pilots, whereby thrust is trimmed to attain the 
desired flight path while elevator control is used to stabilize the aircraft and maintain speed. In addition, 
during operation on the back side of the power curve, this technique is essential because flight path 
response to elevator is unstable for constant thrust while the relationships: 
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The system shown in Figure 1 represents the original conceptual design which incorporated design 
features and characteristics identified during the previous autothrottle/autopilot improvement studies. 
These characteristics are described in the following paragraphs. 

5.2.1 Flight Path Inner Control Loop 

The flight path inner control loop uses flight path error (ye) to derive the throttle command signal as 
flight path angle is directly related to thrust: 

AT = gAy 

Altitude and altitude rate control modes are simply outer control loops which generate a y CMD signal to 
drive the inner loop. 

5.2.2 Inner Speed Control Loop 

The inner speed control loop uses flight path acceleration, which is again related to specific thrust: 
ATocV 

True airspeed control forms the outer loop to the flight path acceleration loop. The flight path 
acceleration command is formed by multiplying airspeed error by the outer loop gain. As with the flight 
path loop, additional airspeed related control modes can simply be added to the front end of the airspeed 
control mode. For example, Mach control can be provided by simply converting the change in Mach 
command to a change in true airspeed. 
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5.2.3 Flight Path Angle Error Signal 


The key design feature is the calculation of the potential flight path angle error signal y £ . This signal 
may be considered as the energy that must be added to the system in order to restore the nominal height 
and airspeed (e.g., in response to a wind gust). The signal is shaped through a high pass filter and used 
for short term control of both the elevator and thrust. In constant energy airplane dynamics (e.g., 
phugoid type motion) the potential flight path angle error signal (y e ) would be zero. 

5.2.4 High and Low Pass Filters 

For large term control, the high and low pass filter yield a unity transfer function for y e , which is then 
integrated to develop the trimmed throttle position. Similarly, in the elevator path, the filters give a 
unity transfer function to V , which is then integrated to develop the trimmed elevator position. 

5.2.5 Pitch Stability 

Short term pitch stability is provided by a conventional 0/0 inner loop. 

5.2.6 Outer Loop Gains 

The outer loop gains of altitude and speed are selected to provide a good stability and transient response. 
Further, for speed control, the ratio of airspeed error gain and inertial acceleration gain is selected to 
achieve maximum cancellation of turbulence induced signal components. 

5.2.7 Cross Coupling Signal Paths 

The cross coupling signal paths are the same for all control modes in order to achieve consistent control 
decoupling for each mode. 

A stability analysis of this system showed that the system was unstable due to the sign of The sign 

must be negative for short term control, but positive for steady state condition. Therefore, the short term 
elevator and throttle command processing require separate y c and V c signal inputs. Crossfeed signals 
summed downstream of integrators would serve for short term control. 
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This concept was developed so that the throttle could be used as an energy controller where: 

Et £ = gY e + V £ 

In addition, it was considered practical to use and y & to develop the short and long term elevator 
command. When the total energy error is zero (i.e., the acceleration error V c is equal to -gy c ), the signal 
inputs and -gy e ) are equivalent and contribute equally to the elevator command. 

This reconfigured system was optimized for high speed (i.e., cruise configuration) as shown in Figure 2. 
The proportional total energy signal to the throttle has been discarded because a stability analysis 
indicated that this signal was not required for adequate stability. 

The stability characteristics were: 


Q_ 

L 

6.7 

.998 

6.16 

.719 

.097 

.888 


5.3 LOW SPEED CONFIGURATION 

The gains derived for the cruise configuration resulted in a slow transient response to command signals 
at low speed. Root locus analysis techniques were used to obtain an improved response. Figure 3 shows 
the resulting configuration and gain values. The primary modifications were: 

(1) Increased system bandwidth by increasing gains Kh and K v . 

(2) Increased inner loop gains for improved stability. 

(3) Change in KDEI, KUD, and KHD gain values to improve cross coupling errors. 
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HCMD~ft h~fps 
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Figure 2. Cruise Configuration 
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Figure 3. Initial Low Speed Control Law Configuration 
















The stability characteristics of the system were: 


L 

C0_ 

.439 

3.91 

.743 

.285 

.969 

.116 


Typical transient responses are shown in Figure 4 for a 100-ft step, Figure 5 for a 10-kn step, and Figure 6 
for a 3° step in flight path angle. 

The double maneuver is of particular interest in examining the performance of the integrated elevator 
and throttle control law. Figure 7 shows the time response of a conventional autopilot (in altitude hold 
mode) and autothrottle for a simultaneous commanded 10 kn increase in airspeed and 105 ft decrease in 
altitude. In this maneuver the throttle moved forward 6° to null the speed error, causing a slow response 
with overshoot. The same maneuver for the integrated system (fig. 3), shown in Figure 8, gave a 
smoother response with no altitude or airspeed overshoot. With this maneuver the throttle settled 
quickly after an excursion of only 2°. 

5.4 WIND SPEED AND TURBULENCE PERFORMANCE 
FOR LOW SPEED CONFIGURATION 

The performance of the baseline integrated autopilot/autothrottle in turbulent conditions is shown in 
Table 1 for 1 fps rms longitudinal and vertical gusts. For comparison purposes, the results obtained for 
the TC V glide slope track law and autothrottle are also presented. It can be seen that the baseline system 
has worse path tracking and elevator activity, although airspeed tracking is better particularly in 
vertical turbulence. The performance of this system in wind shear ( 1 kn/s horizontal shear) is shown in 
Figure 9. The 10:1 ratio of inertial acceleration to airspeed causes a 10 kn steady state airspeed error. 

Three design modifications were evaluated to improve this performance: 

(1) Complementary Inertial/ Airmass referenced acceleration 

(2) Complementary Airmass/Thrust derived acceleration 

(3) Adding a two-dimensional wind stress detector 
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Table 1. RMS Values for Dryden Spectrum Input 


NOISE 

INPUT 

OUTPUT 

BASELINE 

SYSTEM 

V,/V A IR 

X = 20s 

v t ,y 

T = 20s 

v t>y +wsd 

T = 15s 

TCV 

GLIDE SCOPE* 
TRACK + ATT 


V e 

.9831 

1.013 

.9931 

.9876 

1.025 


H 

.6490 

1.065 

1.0726 

1.1826 

.4181 

U 

8 

^THR 

.4184 

.6383 

.6367 

.9380 

.4047 


8 e 

.3662 

.1357 

.3319 

.7410 

.1672 


V e 

.1647 

.1574 

.5485 

.4973 

.4486 


H 

1.9151 

1.973 

.9535 

1.043 

.7916 

w 

8 

i 

^THR 

.8089 

.8082 

.5224 

.6703 

.7795 


5 e 

.7056 

.6966 

.2406 

.2282 

.3116 


* No wind shear detector 


The steady state error due to wind shear can be eliminated by the addition of a complementary filter 
using true airspeed (figs. 10 and 11). Inertial acceleration is filtered through a high pass filter and 
complemented with a lagged airspeed rate. However, this improvement in wind shear performance is 
achieved at the expense of throttle activity (Table 1). 

An alternative method of obtaining acceleration was considered to determine whether throttle activity 
could be reduced. This method consisted of deriving an equivalent acceleration signal from engine thrust 
(EPR) and flight path angle feedback. 


The equivalent acceleration signal was approximated by: 

_g_ AT - gAy 




CO 


Drag effects have been neglected. The thrust can be derived from EPR or engine speed (N x ) and flight 
path angle can be derived by: 

Y - A_ 

vg 

The acceleration signal V T ,y was then used in an identical manner to V G (i.e., complemented with V Air ) 
as shown in Figure 12. 
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Figure 1 1 . Wind Shear Response with Complementary V 

















The overall system stability was decreased by using V T ,y, compared with Vq . The effect of this reduced 
stability can be seen in Figure 13, which shows nearly twice the error of the baseline system (fig. 4) for a 
100-ft step in altitude. Table 2 shows the effect on damping characteristics of varying the V T ,y, system 
complementary filter time constant. 

Table 2. Effect of Varying the Complementary Filter Time 
Constant (X) for V t ,y System 


X = 10 

X = 20 

X = 30 


CO 

<; 

CO 


CO 

.353 

1.82 

.359 

1.84 

.361 

1.84 

.689 

.212 

.724 

.208 

.737 

.207 


The wind shear performance of this system is compared with the performance obtained with the Vj/V AIR 
system in Figure 14. 

Figure 14 shows the effect on rms throttle activity and velocity error due to wind shear of varying the 
complementary filter time constant (x), for both of the systems previously discussed. It can be seen that 
these two systems are almost identical. 

The overall performance of the different control laws is shown in Table 1. The significant difference in 
performance is the effect of vertical turbulence where noise levels on height, throttle and elevator are 
significantly reduced, although velocity error increases by 250%. 

A wind shear detector (WSD) was incorporated with the design (figs. 15 and 16) to complete the analysis 
of wind shear and turbulence. The purpose of this design is to command a new thrust level to offset 
accelerations due to wind. The details of the design are given in Appendix A. 

The rms results (shown in Table 1) for the system with WSD are very similar except for increased throttle 
and elevator activity due to vertical turbulence. However, these results were calculated for a different 
time constant (X). Plotting rms throttle against wind shear for both system shows the characteristic is 
essentially the same, although for a different value of time constant (i.e., the second system with X = 10 
gives similar results to the third system with X =30). (See fig. 17.) 
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RMS THROTTLE POSITION DEG/FPS 
(1 FPS RMS DRYDEN WINDS) 





5.5 INNER LOOP DESIGN 


Though consideration has been given to the performance of this system in wind shear and turbulence, it 
was not considered optimum in terms of stability or path tracking capability. Hence, further root locus 
analysis was carried out to improve stability with the resultant gain configuration presented in Figure 18 
along with the bandwidth and damping of the dominant modes. The performance of this system in 
turbulence was compared with the previous system and found to have lower elevator activity due to the 
lower gain in the proportional h and h paths. However, the revised system gave poorer altitude tracking. 

Previous research work on TSRV auto/and flare control laws had established that K feedback in the inner 
elevator control loop exhibited better path tracking than the conventional 0/0 systems. Therefore, an h/0 
system was derived (see Appendix B). The h/0 system is shown in Figure 19 together with the damping 
and bandwidth of the dominant modes. Altitude and airspeed decoupling were not significantly affected 
by the different inner loop, although wind shear performance was improved slightly (0.7kn for 1 kn/s wind 
shear input). Table 3 shows the effect of turbulence (for 1 fps rms input). The h/0 showed improved height 
tracking compared with the 0/0 loop, although this was achieved at the expense of elevator activity, with 
throttle activity not being significantly affected. 

The performance of the inner loops were compared in the frequency domain (fig. 20a-h). The h/0 system 
showed better low frequency performance for altitude tracking but worse high frequency elevator 
response. In order to take advantage of the preferred characteristics of each system, the two inner lpop 
techniques were combined through a complementary filter. The low frequency component of the h/0 loop 
was combined with the 0/0 high frequency component to achieve the complementation. This was an 
approximate complementation since some of the less significant outer loop terms were omitted to give a 
system consistent with the total energy control law concept. The system diagram and stability 
characteristics are shown in Figure 21. Figures 22 and 23 show the time response of the system to a 100-ft 
altitude command and 10-kn speed command respectively. The altitude and speed decoupling were 
excellent in both cases. In addition, elevator activity had been significantly reduced (e.g., .81° rms to .16° 
rms for vertical turbulence of 1 kn/s rms) without penalizing velocity or path tracking. 

The final inner loop considered for implementation was an angle of attack inner loop. The motivation for 
considering this design was the possibility of providing a simple and effective angle of attack limiting 
capability. 
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Figure 18. Control Law with Improved Stability 
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Figure 21. System with Complementary 9/9, h/6 Inner Loop 
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Figure 22. Time Response to 100-ft Altitude Command 
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Figure 23. Time Response to 10-kn Speed Command 



Table 3. RMS Turbulence Levels (for 1 fps rms Dryden Winds) 


TRANSFER 

INPUT 

FUNCTION 

OUTPUT 

0/0 

(<X AIR) 

h/0 

COMPLEMENTARY 

FILTER 

h/0/0 

a AIR 

U g 

V e 

.9663 

1.013 

.9999 

.9663 

U g 

H 

1.493 

1.063 

1.0262 

1.494 

U g 

8thr 

0.6344 

.633 

.6362 

.6344 

U g 

Selev 

0.1965 

.2143 

.1787 

.1965 

W g 

V e 

0.5038 

.7199 

.6139 

.5610 

w g 

H 

1.3974 

.6587 

.6394 

.9031 

W g 

Sthr 

0.5312 

.5121 

.5174 

.5154 

w g 

Selev 

.0972 

.8084 

.1597 

.7964 


The total energy control system has a proportional flight path angle signal to the elevator, and is easily 

converted to an angle of attack inner loop. The practicality of this depends on the feasibility of accurately 

measuring the angle of attack. The conversion of the oc inner loop can be accomplished by substituting 0 

= y + a for 9 in the inner loop of the 0/0 system, 0 = 57.3 _h + « and for small angles. This system 

V 

is shown in Figure-24. G 

The two formulations for angle of attack (i.e., inertial referenced OC ( OC IN ) and air referenced a ( a A i R )) 
were investigated for the analysis. In the linear analysis, the low frequency signal ( a IN ) was not directly 
affected by longitudinal gust. In addition, the a AIR was linearized in such a manner that it was 
dependent only on vertical gusts. 

The stability characteristics (for both <* IN and a AIR ) are the same as that of original 0/0 system, except 
for the vertical turbulence results for the air system were different to the 0/6 . The <* IN system results 
were identical with the 0/6 system. These results show that a a AIR feedback system improves path 
tracking but at the expense of elevator activity. Elevator activity increased to 0.79° rms for the a AiR 
system compared with 0.097° rms for the system. 

Although the h/0 inner loop offered the best overall performance (i.e., low throttle activity with minimum 
elevator activity) it was decided to proceed with the inner loop. The primary criterion for this decision 
being the possibility of providing a simple, effective °c limiting capability. This concept is reported in the 
Section 6. 
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6.0 DESIGN CONFIGURATION IN NON LINEAR OPERATION 


In a situation which demands a large decrease in velocity, the angle of attack can reach unacceptable 
values and must be limited to prevent stall. Similarly, in response to a command for a large increase in 
flight path angle, the throttle or EPR can limit. In order to cope with these situations, it was necessary to 
modify the linear design. 


6.1 THROTTLE OR EPR LIMITING 

The command cannot be satisfied when the throttle or EPR reaches its limit in response to a large flight 
path angle (y ) command. Localise ol the ci ossleed A (lig. 2J», the error in Light path angle is tied to the 
inner elevator loop causing a bias signal in speed control. Figure 26 shows the result of a 10° flight path 
command from the trimmed condition. In this case, the throttle limited and only achieved about 6°. The 
bias signal (Yqmd " V) was fed to the elevator inner loop and acted as a V cmd signal causing an increasing 
error in airspeed. 

The obvious solution to prevent the problem of bias in speed for throttle limiting was to cut crossfeed A. 
However, simply switching out the crossfeed caused undesirable transients. Therefore a high pass filter 
(washout filter) was added (fig. 27) to eliminate these transients. The action of the filter can be seen in a 
demand in flight path angle (fig. 28). Whenever the throttle limited, switch 1 closed and the crossfeed 
signal tended to zero. In situations when the throttle came off the limit, switch 1 opened and the system 
returned to linear operation without severe transients (fig. 29). 

6.2 ANGLE OF ATTACK LIMITING 

There was a requirement to restrict angle of attack to a suitable small angle to prevent stall. One 
solution was to limit the inner loop command. However, it was found that in response to a large step 
decrease in velocity, the maximum angle of attack limit was reached but the speed command was not 
satisfied. Consequently, a bias was supplied to the engine via crossfeed B (fig. 25) which prevented the 
control of flight path angle. 

The simple solution used in paragraph 6.1 cannot be applied to the case of limiting. Simply switching out 
the crossfeed B resulted in a system where altitude is controlled with thrust. This system was unstable 
with the configuration shown in Figure 25 for realistic values of gain KH. Figure 30 shows the effect on 
airspeed of switching out crossfeed B and holding the system at an a limit of 4°. 
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Figure 25. Tecs Linear Cross feed Model 
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Figure 27. Addition of Washout Filter in Crossfeed A 
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Figure 28. Effect of Washout Filter— Throttle Limiting Due to ^CMD 
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Figure 29. Throttle Limiting Due to Fl cmcj 
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The problem was solved by switching the outer control loop, when limiting occurs, from speed control to a 
control. Figure 31 shows the control configuration. When signal A was greater than B, the system 
switched in a transient free manner and the outer loop became an °c control system. 


K<x _ dV 
K y doc 

Where 

dV Rate of change of V 

doc ~ W.R.T.oc 


This ensured that the velocity and oc loops would produce an equivalent inner loop signal. 

Figures 32 and 33 show the results achieved using the nonlinear system. For a 20-kn decrease in true 
airspeed from 120 kn (fig. 32), the system limits at 4° (fig. 33), and a steady state error of 4 kn remains in 
airspeed. 


6.3 ENGINE CONTROL LOOP 

A control system was developed to satisfy the thrust response requirements of the autothrottle functions. 
This system incorporated a proportional feedback control loop designed to improve the dynamic and 
steady state performance of the 737 engine. In addition, EPR and throttle limiting circuitry was added. 

A simplified block diagram of the control system is shown in Figure 34. 

The JT8D-9 turbofan provided the basic engine for the simulation. The model was adapted from the 737- 
200 model available on the Harris Flight Simulator (Reference 1) and developed into an ACSL 
simulation. 

The input to the engine model was throttle lever angle. The throttle was driven by a servomotor modelled 
by a lag having a time constant (Tg) of 0.15s and rate limit of 10°/s. In addition, hysteresis was added to 
simulate cable backlash in the linkage. The magnitude of the hysteresis was obtained from Reference 2. 

The proportional feedback control loop was designed using a linear engine model consisting of a gain and 
a simple lag Te = Is. Figure 35 shows that the control loop consisted of a second order system comprising 
of engine model and throttle servo. Assuming that the steady state gain should be unity and a damping 
ratio (£) of 1, then values were calculated for KEPRP (1 12.5) and KFB (.55). 
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Figure 31. System Configuration Induding a Outer Loop 
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Figure 34. Engine Control System 

















EPR Demand Engine Loop (Linear Design) 






The simulation results obtained using these gain values are shown in Figure 36 for the initial conditions 
of: 

Alt = 20000 ft 

Mach = 0.69 

EPR, C = 1.55 

The commanded change in EPR is 0.5. Owing to backlash (of 2.17), the actual EPR started below the 
commanded level. It can be seen that the rate limit within the throttle of 107s greatly modified the input 
command to the engine (TAUTO) when compared with throttle command (THR CMD). The purpose of 
the engine control loop is to provide the required thrust in response to the outer loop command signal. 
However, it is not possible to feedback thrust measurements and EPR feedback (normally employed with 
Pratt and Whitney engines). A linear relationship between EPR and thrust exists provided atmospheric 
conditions remain constant. Nevertheless, variation in altitude will change this relationship. In this 
system an attempt has been made to correct this change by modifying the outer loop command via 8 (ratio 
of P/Po). Furthermore, during flight the aircraft weight will decrease due to fuel consumption modifying 
the required thrust command. Figure 34 shows a time varying gain (W/Wo) has been introduced to 
compensate the EPR command for weight variation. 

Two important features of the overall engine control systems are the throttle limit and EPR limit 
feedbacks (fig. 37). If either the EPR command or throttle command exceed the 737 engine limits, then a 
scaled signal is fed back to reduce the integrator input and prevent further saturation of the integrator. 
The action of the limits is demonstrated in Figure 38 for a step change in engine outerloop command (not 
shown in Figure 38). The output of the integrator (ICTI) ramps up until the throttle limit is reached (65°), 
at which point a signal (DELTA) is fed back to reduce the integrator input. It can be seen that the throttle 
command (THCM) and throttle output (TAUTO) are held at the 65° limit. 
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Figure 37. Detailed Engine Loop (ASCL Simulation) 


























7.0 ADDITIONAL CONTROL MODES 


The integrated autopilot/autothrottle design, discussed in previous sections, was configured to control 
height (h), true airspeed (V TAg ) and flight path angle (FPA, Y). However, the basic design philosophy was 
such that the inner loops are driven by FPA error (y c ) and longitudinal acceleration ( V ). It is therefore 
relatively simple to incorporate additional control modes by designing these modes to generate the 
appropriate y e and V signals. 


7.1 GROUND SPEED MODE 

Control of ground speed was the simplest additional mode to implement. Implementation consisted of 
providing a ground speed command, and feedback of ground speed V c . and ground speed rate V ( . in place 
of the true air speed signals (fig. 39). 


7.2 MACH AND CAS SPEED MODES 

The Mach and CAS speed control modes are required to facilitate profile descent procedures. One 
preferred procedure is to descend at constant Mach to some preselected altitude at which point switch to 
descend at constant CAS. Conversely, during ascent one procedure is to fly at constant CAS then switch 
to constant Mach. 

The means of implementing the Mach/CAS control loop (fig. 40) is documented in Reference 3. This 
implementation consists of converting the MACH or CAS signal to a TAS signal which drives the velocity 
control loop. 

In the case of the Mach control law, the Mach command signal is converted to V TCMD by computing the 
current speed of sound and using this as the conversion factor. In the case of CAS, the V CAg error signal 
is formed prior to converting to V TAg . This technique allows an approximate conversion factor (derived 
from speed conversion tables) to be used (i.e., V^ Ag = V 1 (1 - .12 x 10 4 h). 

It can be seen that as V TC is controlled to zero V CAg£ must go to zero, and that an exact measurement of 
height is not critical to accurately controlling V CAg . During the first portion of descent, the control 
system is designed to hold the airplane at constant Mach. When CAS equals a preselected CAS, then the 
system automatically switches to CAS hold and continues the descent by maintaining constant LAS. On 
ascent the reverse is true. 
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Figure 40. Implementation of Mach/CAS Hold 













7.3 GLIDE SLOPE AND VERTICAL PATH MODES 


An additional control loop has been added to the basic configuration to allow a transient free captive 
exponential of the glide slope during landing approach. Figure 41 shows the geometry of the engagement. 
The airplane is approaching from beneath the glide slope at a current height error h e (ft) off the glide 
slope. An exponential capture law dictates that the airplane follow a trajectory defined by 

^derived 

where 

^derived = ^ res P ect to glide slope 

This law was accomplished by forming an h derived signal by complementing the low frequency component 
of h e rate with the high frequency component of hj using the appropriate filtering (fig. 42). For an h e that 
is positive and decreasing (i.e., the airplane is approaching beneath the glide slope), the error rate h derived 
measured with respect to the glide slope is positive (fig. 42). The effect of this signal is that the output of 
glide slope mode (y £ ) is 'biased' such that it changes sign and commands a negative flight path angle prior 
to the glide slope. Obviously the distance h e at which V change sign is a function of h derived and forward 
gain (for constant airspeed) as given by the previous equation: 

h = derived 

K 

D X ^derived 

k h 

Therefore engagement of glide slope mode will occur earlier for high h derived or low K h . 

Vertical path mode has been designed to function in a similar manner to glide slope (fig. 43), however, the 
control input is generated by the flight control computer in the form of a desired height profile. 

Additional software has been added so that several vertical paths can be flown consecutively with 
transient free switching being achieved between each leg. 
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Figure 42. Transient Free Switching— Glide Slope Mode 


















8.0 MODE CONTROL PANEL 


In addition to developing the integrated autopilot/autothrottle algorithm, consideration has been given to 
the hardware necessary for a pilot to select the various autopilot modes and the software necessary to 
control the arming, engagement and switching between the various modes of operation. Figure 44 shows 
the mode control panel designed for the system. 

Figure 44 shows that the panel is split into four sections: longitudinal modes, lateral modes, speed modes 
and system status. At present the lateral modes (i.e., longitudinal V-CWS) and PROFL (speed profile) 
modes have not been designed. However, both are intended to be included at a future date. Pushing the 
appropriate mode control button changes or arms the modes. A longitudinal mode must be selected in 
addition to a speed mode. Engagement of a mode causes the corresponding mode button to light up green, 
whereas arming a mode causes automatic engagement to occur once certain criteria have been satisfied. 
In this case the mode control button lights up amber when armed and turns green when engagement 
occurs. The operation of the individual control modes is discussed in the following paragraphs. 

8.1 LONGITUDINAL MODES 


8.1.1 Flight Path Angle 

The flight path angle control mode is engaged without prior arming. (Engagement causes tracking of 
commanded FPA.) When disengaged, the FPA display shows actual FPA. However, when engaged the 
actual FPA is taken to he the command and the display. Once engaged, the command can he adjusted 
using the adjust knob. 


8.1.2 Velocity Control Wheel Steering (CWS) 

Velocity CWS is a computer augmented manual control mode in which the rate of change of commanded 
FPA is proportional to control column deflection. The mode is engaged directly and is compatible with 
any of the speed modes. 

8.1.3 Glide Slope (GS) and Vertical Path (V-Path) 

GS mode is armed by pushing the GS mode button with engagement occurring automatically as described 
in paragraph 7.3. The mode will not arm without detection of a valid GS signal. V-PATH mode operates 
in a similar manner to GS. However, in this case, the mode will not arm unless a vertical path is stored in 
the flight control computer. 
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8.1.4 Altitude 


The altitude hold mode can either be engaged directly or armed by preselecting a desired altitude. When 
the mode is not engaged or armed the display shows the current altitude in feet. Engagement causes the 
actual altitude to be fixed on the display and to be taken as the command signal. The display can be 
altered using the altitude adjust knob. 

The adjust knob can also be used to establish the altitude preselect mode by dialing the adjust knob when 
the mode is disengaged. This causes the existing altitude at the time of knob rotation to be stored in 
memory. The altitude increments are added to the stored altitude to hold the new preselected altitude. 
The mode is then armed and will automatically engage when the capture criterion is satisfied. 

8.1.5 EPR 

The EPR mode causes the system to command maximum safe thrust setting while holding speed. This 
mode can be engaged directly. Selection of another flight path mode while the EPR mode is engaged 
causes the EPR mode to disengage and the selected FP mode to engage. An additional pushing of the EPR 
button following engagement, will cause the EPR mode to disengage and the FP mode to engage. 

8.2 SPEED MODES 

The speed control options on the mode control panel (fig. 44) are: CAS, MACH and PROFL. A speed mode 
will always operate in conjunction with a longitudinal mode with the default mode being CAS. 

Switchover between CAS and MACH is designed to occur automatically at certain flight conditions, but 
may be overridden by pushing the desired mode switch. 

The display shows either the commanded CAS or MACH number (depending on mode selected) which is 
equal to the actual speed at the time of engagement. 

The command can be changed by rotating the adjustment knobs. Changing between CAS and MACH 
causes the equivalent speed command to be displayed, preventing transients on switchover. 

The speed profile mode (PROFL) can be armed by pushing the PROFL button. Engagement occurs 
automatically once the desired criteria are satisfied. For the mode to arm, a speed profile must be stored 
in the flight control computer. 
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9.0 PERFORMANCE ASSESSMENT 


The control law algorithm was implemented on the Harris non-linear airplane simulator to demonstrate 
the overall system concept and evaluate the performance using a more complex non-linear simulation. 
The implemented system is shown in Figures 45-48 with Figure 45 showing a simplified block diagram of 
the total system and Figures 46-48 showing the detailed system. 

9.1 ENGINE CONTROL LOOP 


Comparison runs were made with the Harris simulator to confirm that the ACSL simulation (paragraph 
6.3), had yielded a realistic engine response. These were made at two aerodynamic conditions (figs. 49 
and 50): (1) altitude = 5000 ft, equivalent air speed = 150 kn, and (2) altitude = 20000 ft, equivalent air 
speed = 310 kn. 


Gains predicted by linear analysis resulted in an underdamped response with excessive throttle 
excursion, and consequently, the forward and feedback gains (fig. 46) were reduced to KEPRP = 75 and 
KEFB = 0.3. The effect of reducing the system gain did not significantly effect the response time of the 
engine to the step input ( DEPRO = 0.5). The command signal (DEPRO, fig. 49) is defined as: 


EPR Cmd - EPR min 

where EPR . = 0.99 

min 


Hence, the output EPR changes from 1.4 to 1.9 for a command from 0.4 to 0.9. Figure 49 also shows the 
throttle command (THCM) and throttle (TAUTO). It can be seen that a dominant feature of the transient 
response is the rate limiting of the throttle (10%). Figure 50 shows a 0.5 change in EPR from an initial 
value of 1.55 at the cruise condition. 


9.2 SPEED AND ALTITUDE MODES 

The primary operating mode of the system is in altitude and speed hold modes. To demonstrate the 
system performance with these modes engaged, the response to speed and altitude steps was recorded 
over a range of aerodynamic conditions: 

(1) 120 kn CAS, 1500 ft (figs. 51 and 52) 

(2) 150 kn CAS, 5000 ft (figs. 53 - 55) 

(3) 200 kn CAS, 10000 ft (figs. 56 and 57) 

(4) 250 kn CAS, 15000 ft (figs. 58 and 59) 

(5) 320 kn CAS, 20000 ft (figs. 60 - 62). 
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Figure 45. Integrated Elevator/Thrust Control System 
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Figure 46. Engine Control Loop 
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Figure 48. Thrust and Elevator Command Processing 
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Figure 49. Engine Performance (5000 ft, 150 kn, EAS) 
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Figure 59. 100-ft Step in Height (250 kn, 15000 ft) 
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Figure 60. 10-knStep in V QAS (320 kn, 20000 ft) 
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Figure 61. 100-ft Step in Height (320 kn, 20000 ft) 
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Figure 62. Double Maneuver ( AVcas = 10 kn, A H- -380 ft) 


Height error was less than the design requirement of 20 ft for all the velocity changes shown. A typical 
height error for 10 kn change in velocity was approximately 10 ft (fig. 53). The transient response in all 
cases exhibited no overshoot and reached 95% of final value by about 35 s. In all cases the engine 
response was smooth with no excessive or unnecessary throttle motions. Vertical accelerations peaked at 
O.lgfor 10 kn changes in V CAS , except at condition (4) (250 kn CAS, 15000 ft) when the 'g r level reached 
0.15g. 

For all height changes shown, velocity error did not exceed the design requirement of 1 kn and was 
typically 0.5 kn. As with velocity changes, the transient response exhibited zero overshoot. 

In Section 2 is discussed one design feature of the integrated autopilot/autothrottle (i.e., the constant 
energy concept). In a constant energy maneuver (i.e., one in which the change in potential energy is 
matched by the change in kinetic energy) the system will use the elevator to retrim the energy in the 
system and the throttle activity will be negligible. This type of maneuver is demonstrated in Figure 55 
where the command inputs are + 10 kn change in velocity coupled with a -138-ft change in height. These 
command inputs were calculated to be energy equivalent at 150 kn. It can be seen that the elevator 
immediately responds at the start of the maneuver, whereas throttle activity is negligible. The throttle 
command (THCM) does show about 2° change over the maneuver, but owing to backlash in the system the 
throttle lever angle (TAUTO) does not respond with no change in EPR detectable. 

A similar energy equivalent maneuver was performed at 20000 ft with a 10-kn change in speed 
corresponding to a -380-ft change in height. Again the change in throttle command is small ( < 4°) and 
the change in EPR and throttle lever angle negligible. 

The effects of large changes in command inputs are shown in Figures 63 and 64. Figure 63 shows the 
response to a 1000-ft change in altitude and demonstrates the action of two of the limiters in the system. 
These limiters are designed to limit the rate of climb and the acceleration levels experienced during large 
maneuvers. 

The reference to Figure 47 shows a limiter in the altitude command path. The system is designed such 
that a command for change in altitude generates a flight path angle error (y c ) signal which is limited. 
The result is a height rate limit which varies as a function of velocity from approximately 1800 ft/min at 
low speed (120 kn) to 6800 ft/min at cruise (310 kn EAS). 

In addition to this limit, other limiters exist on y c and V £ in order to restrict acceleration commands to 
approximately O.lg. However, during transients, the normal acceleration will exceed this O.lg level. 
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Figure 64. Large Velocity Change (20 kn change 150 kn, 5000 ft) 
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Figure 63 shows that during the 1000-ft altitude change, the height rate was limited to approximately 30 
fVs. During the initial phase of the maneuver the normal acceleration peaked at 6 ft/s 2 but settled to 4.4 
ft/s 2 . Capture of the new altitude was achieved with no overshoot and velocity error was insignificant 
(< lkn). 

The effect of a 20-kn step in CAS (fig. 64) shows negligible height error and smooth velocity captive. 
Transient peaks in normal acceleration exceeded O.lg reaching 6 ft/s 2 . However, future adjustment of 
limiters or gains could reduce this figure if necessary. 

The speed mode options available with the integrated autopilot/autothrottle are CAS, Mach or ground 
speed modes. The design and operation of these modes is described in paragraphs 7. 1 and 7.2. 

Considering the CAS/Mach mode operation, a switch over between the modes was designed to occur 
automatically at preselected speeds. During ascent, the system will switch from CAS to Mach mode 
whereas during descent the system switches from Mach to CAS. 

Figure 65 shows the system operation: initially the airplane is commanded to descend at a flight path 
angle of -3° but holding Mach 0.7. CAS increases to the preselected switch point of 325 kn, then the 
system holds CAS constant while Mach decreases. An added feature of the altitude mode is the capability 
to preselect desired altitude. The system climbs or descends at the required flight path angle until the 
switching criterion is satisfied then engages altitude mode in a smooth transient free manner (fig. 66). 

9,3 FLIGHT PATH ANGLE, GLIDE SLOPE AND VERTICAL PATH MODES 

The response to a step change in Flight Path Angle (Ay = 3°) is shown in Figure 67 (150 kn EAS, 5000 ft) 
and Figure 68(310 kn EAS, 20000 ft). Captive was smooth with minimum change in throttle. However, 
normal acceleration peaked at 8 ft/s 2 in the high speed example. 

The design of glide slope and vertical path following mode are described in detail in paragraph 7.3. The 
aim of the mode is to provide a smooth transient free captive of either the glide slope or some portion of 
the precalculated vertical path. The operation of both these modes is illustrated in Figures 69 and 70 
showing the desired transient response of y and height with low acceleration (AZCG) and negligible 
velocity error. 


9.4 NONLINEAR OPERATION DUE TO « OR ENGINE LIMITING 

The integrated autopilot/autothrottle system has been designed to prevent stall due to excessive <* 
command and to prevent engine overboost due to large throttle commands. The design and 
implementation of these features is described in Section 6. Figure 71 shows the effect of commanding a 
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large flight path angle (y cmd - 10°) when the maximum throttle angle has been limited to 40°. In this 
case, the actual FPA settles at 5° and the elevator retrims to maintain velocity. Maximum velocity error 
was 1.6 kn in this limiting condition. 

The action of the limit is demonstrated in Figure 72a and b. In Figure 72a a 20-kt change in V g is 
commanded causing the angle of attack («) to reach in excess of 1 1°. The FPA deviates by a maximum of 
approximately 1.2". 

The result of setting « f (i.e., maximum allowable «) to 9° is shown in Figure 72b. FPA hold is 
maintained with approximately 1° error, but velocity is prevented from reaching the command input as « 
is limited at 9°. 

The integrated autopilot/autothrottle includes an EPR max mode. This option is designed as a safety 
mode to use as a "Go Around" feature. Engagement of EPR max mode causes the engine to deliver 
maximum thrust as quickly as possible. Velocity hold mode is still engaged, so the airplane climbs at 
maximum height rate with the limits on height rate and acceleration not applying. Figure 73 shows the 
EPR max mode in which the throttle ramps forward at maximum rate with EPR limits at the maximum 
safe value (2.25). Elevator activity is considered high and velocity error reaches 5 kn. However, this is an 
emergency mode and these factors were acceptable. 

9.5 TURBULENCE AND WIND SHEAR 

The effect of wind shear and turbulence on this control law was investigated using the Harris simulator. 
The Velocity Control Loop employs a complementary filter to provide the derivative feedback signal (see 
paragraph 5.4). Simulation runs were made with horizontal turbulence (Dryden wind spectrum lft/s 
rms) to examine the result of variation of the complementary filter bandwidth (CO R = l/ tk ) . 

The standard deviation (SD) of several system states are shown plotted against C0 k in Figures 74 and 77 
for two aerodynamic conditions: (1) H = 20000 ft, V = 310 kn, and (2) H = 1500 ft, V = 120 kn. 

It is apparent that the SD of elevator (og e ) in Figure 74 and throttle (° THCM and ° TAUT0 ) in Figure 75 
were very sensitive to variation in C0 K . An increase in C0 K caused an increase in the SD. The SD in 
throttle command (° THCM ) was larger than o TAUTO (throttle lever angle) due to backlash. The SD of 
height was sensitive to variation in CO K at high altitude and speed, but remained at 2 ft ± 10% at the low 
speed and altitude. Variation in velocity SD was small. In vertical turbulence (CO g = lft/s rms) the 
response was largely independent of variation in co^ (Table 4) at both aerodynamic conditions. 
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Figure 75. Horizontal Turbulence —SD of Throttle Against 
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Table 4. Effect of Vertical Turbulence (1 fps rms) 
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The performance of the system was examined for horizontal wind shear (modelled as a 1 kn/s ramp). 
Maximum velocity and height error were plotted for various valves of CO K at the two aerodynamic 
conditions (high and low) previously described. Figures 78 shows that increasing the filter bandwidth C0 K 
( i.e. , increasing the bandwidth of the velocity loop) reduces the maximum velocity error. 

However this is achieved at the expense of height error which increases significantly with increase in 
bandwidth at high altitude (fig. 79). The effect was far less marked at low altitude and speed. An 
interesting comparison can be made between the wind shear and turbulence performance: Maximum 
error and SD of height increase with increasing filter bandwidth, whereas maximum velocity error vary 
significantly with filter bandwidth, with the SD of velocity insensitive at high altitude and only 
decreased slightly at low altitude. 

The linear results plotting SD of throttle against maximum velocity error V for varying filter bandwidth 
(paragraph 5.4) shows that a trade-off is possible between throttle activity and velocity error wind shear. 
Figure 80 shows that throttle lever angle (TAUTO) activity was lower both at high and low speed than 
the linear predictions. However, the throttle command signal was approximately 75% higher than the 
linear prediction. To take advantage of the ability to trade off throttle activity in turbulence against 
maximum velocity error in wind shear, the system time constant (X K ) has been programmed to vary as a 
function of altitude. Thus low altitude (1000 ft and below) X K = 4, thereby reducing the velocity error 
due to wind shear. At higher altitude the effect of wind shear is not critical and it is desirable to minimize 
throttle activity. Therefore above 2500 ft X K is fixed at 20. Between 1000 to 2500 ft, X K varies in a linear 
manner with altitude. Details of the implementation are given in Figure 47. 

The turbulence and wind shear results did not show a significant difference when either V G or the 
^derived s ^S na l was implemented on the Harris simulator. Further efforts in modelling the 
instrumentation required to determine thrust and current weight are needed before an accurate 
evaluation of the two methods can be made. 
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Figure 78. Horizontal Wind Shear— Velocity (Vq) Against Wk 
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10.0 CONCLUSIONS 


1. The integrated autopilot/autothrottle was designed to give: 

(a) Fuel efficiency by employing the total energy control concept and minimizing throttle 
activity. 

(b) Low development, implementation and certification costs by designing the system around 
a fixed inner loop configuration. 

(c) Maximum safety by preventing stall and preventing engine overboost. 

2. The original design constraints have been satisfied during normal operation (e.g., maximum 
velocity error during altitude maneuvers was less than 1 kn and maximum height error during 
velocity changes was less than 20 ft over the full aerodynamic range). 

3. Wind shear performance was improved to satisfy the design requirement of 5 kn maximum error 
with the addition of a complementary inertial/air mass referenced acceleration. However, this 
improvement was achieved at the expense of throttle activity. 

4. Variation of the complementary filter time constant, enabling the system to be optimized for 
wind shear criterion, was satisfied during approach and landing, but throttle activity was 
minimized during cruise. 

5. The use of a 2D wind shear detector did not improve wind shear performance when compared 
with the preferred system. 

6. Several inner loop designs were evaluated. Of these designs, linear analysis showed the 
optimum was an H/G/q complementary inner loop in terms of decoupling in altitude and speed, 
low elevator activity, and best path velocity tracking in turbulence . 

7. A feedback inner loop was selected in an attempt to provide a simple method of limiting. 
However, this method was not satisfactory for limiting, and stall prevention was incorporated by 
an additional outer loop control. 

8. Future work should consider reimplementing the fi/0/0 inner loop. 
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9. An engine controller was designed to satisfy the thrust response requirements of the 
autothrottle functions. Overboost of the engine was prevented by limiting the EPR command 
signal. In addition, the engine controller has an EPR MAX mode designed as an emergency 'Go 
Around’ feature. 

10. Additional modes have been added to control ground speed (i.e., Mach or CAS) and capture the 
glide slope or follow a preprogrammed vertical path. 

11. A mode control panel was designed to enable mode selection by the pilot. Consideration was 
given to the software design necessary to control arming, engagement and switching between 
the various modes. 

12. The overall system concept and performance were evaluated using a more complex nonlinear 
simulation implemented on the Harris airplane simulator. 

13. The Harris simulator confirmed the results achieved by linear design. 

14. Future work should compare the design concept and results of this conventional design with 
results obtained by solving optimal control problem using a suitable cost function. This 
approach may yield further insight into minimizing the effect of wind shear and turbulence. 
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APPENDIX A. TWO-DIMENSIONAL WIND SHEAR DETECTOR 


The basis of the wind shear detector is to calculate a new thrust level to compensate for accelerations due 
to the wind shear. 

It can be shown that the vertical wind is: 


Vw v = -h + V T tan 0 - V T tan « cos (J) (A. 1) 

where h — inertial altitude rate 
V T = true air speed 
0 = pitch altitude 

X = true angle of attack 
(J) = roll altitude 

In addition it can be shown that the change in flight path angle due to vertical wind, assuming constant 
speed and power setting is: 

Ayw v = ~Vw v /V T (A. 2) 


Hence Ayw v = Yi + tan 0 - tan « cos (J) (A. 3) 

or for small angles: 

Ayw v = Yi + 0- a (A.4) 

The thrust required to compensate for such a change in flight path angle is: 

ATw v = WAyw v (A,5) 


Where W = Weight (lbs) 

Equation 4 is implemented by measuring longitudinal acceleration (a ) using a body mounted 
accelerometer, where: 


a x = V I + g 0 


(A.6) 


where Vj = velocity rate (inertial) 
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Figure A.l shows that a lagged 0(0) can be derived by integrating the accelerometer measurement, 
feeding back V TAS , and high pass filtering the resultant signal. This concept is used in Figure A. 2 to 
obtain the AyW v signal. This signal (AyW y ) is rate limited to filter the high frequency turbulence. 
Finally the signal input to the throttle was divided by the throttle-to-thrust gain to get the command in 
terms of degrees of throttle. 

It should be noted that the computation of AT W is dependent on the accuracy of the angle of attack signal 
and that accurate determination of a has not been fully investigated. 
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APPENDIX B. DERIVATION OF h /0 SYSTEM 


Pitch rate may be expressed as: 


0 = — h • 

V Q + a 0, « in rad/sec 

Substitute for the state equation from the aircraft equation of motion: 


(B.l) 


* = (-.3133)U - (.6543) « + (l.)q - (,000723)6 e - (.108x10 ^ 


(B.2) 


Considering small perturbations only from an initial steady state condition, constant thrust, no 

1 

horizontal or vertical turbulence, and the elimination of « with = 0 h, the equation reduces to: 

V G 


1 .. 


1 


9 = 77 ~ h- 33.13U - .6543(0- _ h) + l.Oq - .0007238 

V G V G 


(B.3) 


With 0 = q, the equation can be solved for 0 as follows: 


0 = ,4319h + .2827 h - .0638 o - ,1354V 

e 


(B.4) 


Where V = V G0 U, 0 is in degrees, fi is in ft/sec 2 , 6 g is in degrees and both h and V are in ft/sec. 


From Figure 18, the equation for the elevator command is: 

.75 


8 ec = 4q + 8^0 + g/V G h-l.V 


eng + s 


- 1V TAS- V en e -g/ V G<- lh - h > 


(B.5) 
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Substitution of equation 4 in equation 5 results in: 


8 PC = 4q + 8 


,4319h - ,135V + ,283fi - ,0638 e + g/V G h 


.75 

- V n g + X 


- 1 V TAS- V eng-g /V G(-l h - h ) 


(B.6) 


The elevator position is related to S ec by the elevator servo transfer function: 


Sec = 


ec 


.05S+1 


(B.7) 


Substition of equation 6 into equation 7 and solving for 8 yields: 


ec =- 


.655 


.033S+1 


4q + 8 ^.4319h - .135V + ,283h + g/V G h 
.75 


-V + -Q- 

eng O 


- 1 V TAS- V eng-g /V G(- lh - h ) 


(B.8) 


This equation represents an h/9 control law that has nearly the same stability characteristics as the 
original system. For implementation of this system, the .03 sec lag must be replaced with the actual 0.5 
sec elevator servo transfer function. This would have a very small effect on the aircraft dynamics. 
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